BLOCK-DIAGONALIZATION OF OPERATORS WITH GAPS, 
WITH APPLICATIONS TO DIRAC OPERATORS 



JEAN-CLAUDE CUENIN 



Abstract. We present new results on the block-diagonalization of Dirac 
operators on three-dimensional Euclidean space with unbounded potentials. 
Classes of admissible potentials include electromagnetic potentials with strong 
Coulomb singularities and more general matrix-valued potentials, even non- 
self-adjoint ones. For the Coulomb potential, we achieve an exact diago- 
nalization up to nuclear charge Z = 124 and prove the convergence of the 
Douglas-Kroll-Hefj approximation up to Z = 62, thus improving the upper 
bounds Z = 93 and Z = 51, respectively, by H. Sicdentop and E. Stockmcyer 
considerably. These results follow from abstract theorems on perturbations 
of spectral subspaces of operators with gaps, which are based on a method 
of H. Langer and C. Tretter and are also of independent interest. 
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1. Introduction 

The Dirac operator on R 3 governing the motion of a relativistic particle of 
half-integer spin in the presence of an external electromagnetic field is given by 
(in units where the reduced Planck constant, the velocity of light and the particle 
mass are equal to one) 

(L1) H - V-io-.(V-iA) -1 + $ J' 

Here, er := (ci, 02, 03) is a formal vector whose components are the Pauli-matrices 




The author gratefully acknowledges the support of Schweizer Nationalfonds, SNF, Grant No. 
200020.130184. 



2 



JEAN-CLAUDE CUENIN 



$ : R 3 — > R is the electric potential and A : R 3 — > R 3 is the magnetic vector 
potential, which determine the electric and magnetic field (uniquely up to a choice 
of gauge) by virtue of 

E = V$, B = curlA. 

We regard H as an unbounded operator in the Hilbert space of square-integrable 
functions 

(1.2) H = L 2 (R 3 ,C 2 ) 8 L 2 (R 3 ,C 2 ) . 

It is well known that the Dirac operator is not bounded from below, giving rise 
to an infinite "sea" of unphysical negative energy states. In the context of Dirac's 
hole theory, the Pauli exclusion principle is invoked in order to "fill the sea" , i.e. 
declare the states of negative energy as already occupied and restrict H to its pos- 
itive spectral subspace. While hole theory turned out to be unsatisfactory from a 
theoretical point of view and was ultimately made obsolete by quantum electro- 
dynamics, it still has its merits in atomic physics and quantum chemistry, where 
the energy scale is well below the threshold for particle creation and annihilation. 

In the field-free case (<& = 0, A = 0) the now famous Foldy-Wouthuysen trans- 
formation [15] may be invoked to decouple the positive and negative spectral sub- 
spaces. However, the original method proposed by Foldy and Wouthuysen in the 
presence of external fields makes use of an ill-defined expansion in the inverse 
speed of light and has to be discarded, see [39]. For Dirac operators in a purely 
magnetic field ($ = 0), supersymmetric methods have been employed by Thaller 
[381 139] to construct an exact transformation. However, these methods do not 
apply to electric potentials, which arc of paramount importance from the physical 
point of view. In [29], Langer and Tretter developed an abstract method for di- 
agonalization of block operator matrices based on indefinite inner product spaces, 
which yields an exact transformation of the Dirac operator for bounded electric 
potentials of norm less than one. Siedentop and Stockmeyer [37] proved the exis- 
tence of an exact transformation for the Dirac operator with Coulomb potential, 
as well as the convergence of an approximate block-diagonalization, known as the 
Douglas-Kroll-HcB (DKH) method. The latter, proposed by Douglas and Kroll 
in their seminal paper |14j . consists of a an iterative scheme which decouples the 
positive and negative spectral subspaces up to any given order in the coupling 
constant of the potential. Its usefulness for quantum chemical implementations 
was first realized by B.A. Hcfi in |19| and has since turned into one of the most suc- 
cessful computational tools in relativistic quantum chemistry [20] [33] [34] 135] [36] . 
From a mathematical point of view, the method of |37j was extended to the multi- 
particle case in [52] and to complex-dilated Dirac operators in [3T]. Moreover, it 
was mentioned in |25j that the same method works for non-vanishing magnetic 
fields (A / 0). Different techniques were employed in [17] and in [7] [9] [10] to 
handle weaker electric potentials than the Coulomb potential. 

By using a combination of the methods of [29] and [37] we achieve generaliza- 
tions of the aforementioned results; the main novelties are the following: 



BLOCK-DIAGONALIZATION OF OPERATORS WITH GAPS 



3 



• When $ = —7/I • I is the Coulomb potential and A = 0, we obtain an 
exact block-diagonalization of H up to nuclear charge Z = 124, extending 
[371 Theorem 1], where a transformation was shown to exist up to Z = 93; 

• We show the convergence of the DKH approximation up to Z = 62, ex- 
tending [37l Theorem 2], where the convergence was proved up to Z = 51; 

• Potentials are allowed to be (not necessarily symmetric) scsquilincar forms; 

• The method can be adapted to handle strong (e.g. constant, but also 
unbounded) magnetic fields. 

Moreover, our transformation can be chosen as the direct rotation (see [T2"] ) 
between the subspace of upper (lower) component Dirac spinors and the positive 
(negative) spectral subspace of H, that is, it has minimal deviation from the 
identity among all such transformations. This is a consequence of the fact that 
the positive and negative spectral subspaces Q±H of H admit a representation 
in terms of so-called angular operators. Identifying the direct summands in the 
decomposition Q 1.2)1 with the subspaces 

P U H := I Q : u G L 2 (R 3 , C 2 ) } , := j : g G L 2 (R 3 , C 2 ) j 

of "upper-" and "lower-component" Dirac spinors, respectively, this means that 

Q+n = j : U G P U U\ , Q-H = I ( X ~ V ^ : v G Pi^j , 

where X± are bounded operators in L 2 (R 3 , C 2 ) . In particular, we obtain bounds 
on the norms of X±, which means that for any eigenfunction ip = corre- 
sponding e.g. to a positive eigenvalue of H, we must have 

ll*/||<l|A + ||||* u ||. 

We emphasize that our technique is purely operator-theoretic in nature and 
thus not limited to the Dirac operator. Our main results, Theorems 13.11 13.31 and 
Theorems 13. 71 l3~8l apply to arbitrary self-adjoint operators with a spectral gap at 
zero, perturbed in the quadratic form sense. 

We briefly sketch the outline of the paper. In Section 2, we set the necessary 
operator-theoretic background. Among the tools needed are indefinite quadratic 
form methods, spectral projections for non-selfadjoint unbounded operators and 
a theorem on accretive operators in spaces with an indefinite inner product. The 
main abstract results of the paper are stated in Section 3 and are applied to the 
Dirac operator on R 3 with Coulomb-type potentials (with and without magnetic 
field) in Section 4. Section 5 contains the proofs of the main theorems. 

2. Preliminaries 

2.1. Notation. Let X and Y be Banach spaces. For an operator S(X — > Y), 
we denote by T>(S) C X its domain and by ran(S) C Y its range. All Banach 
spaces are always assumed to complex and operators between them are assumed 
to be linear. The Banach space of bounded operators from X to Y is denoted by 
L(X, Y); if X = Y, we simply write L(X) := L(X, X). For the identity operator in 
X we write Ix or I if it is clear from the context which space is meant. If Y = X, 
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then p(S) and <j(S) denote the resolvent set and spectrum of S, respectively. 
If S is closed, then the former coincides with the set of all z G C such that 
S — z : T)(S) — > X is bijcctive; here, we used the abbreviation S — z := S — zl. 
If S is closable, we denote its closure by S. By an isomorphism between two 
Banach spaces we mean a linear homeomorhism. A subspace C C X is always 
understood to be closed. The topological direct sum of X and Y is denoted by 
X + Y . For two Hilbert spaces H and IC, the orthogonal sum is denoted by H © K,. 
Moreover, the scalar product (•, •) in a Hilbert space H is assumed to be linear 
in the first variable, and for a densely defined operator H in T-L, its Hilbert space 
adjoint is denoted by H* . The Schatten-von Neumann ideals in Ji of order p are 
denoted by S P (H) For a sesquilinear form t>:2?(t>)x2?(t>)— >-<D, the corresponding 
quadratic form is abbreviated by v[u] := v[u,u\. We say that an interval (a, j3) 

is a spectral gap for a self-adjoint operator H if (a,/3) C p{H). An integral J is 
always understood in the sense of the Cauchy principal value at zero and infinity. 

2.2. Indefinite quadratic forms. An operator H in a Hilbert space T-L is said 
to be associated with a densely defined sesquilinear form f) if the following hold: 

i) H is closed and densely defined; 

ii) V(H),V{H*) C V{% 

iii) (Hu,v) = t)[u,v], u e w G V{% 

iv) (u,H*v) = t)[u,v], u G 0(f)), u G Z>(iF). 

If such an operator exists, then it is uniquely determined, see [43l Proposition 2.3]. 
If f) is symmetric and H is self-adjoint, then () is said to be represented by H if 
D(fj) ^Odi?! 1 / 2 ) and 

f)[u,v] = ('|ff| 1 ' /2 u,sign(lf)| J ff| 1 ' /2 i;) , «,« G O(f)). 

The classical first representation theorem as found e.g. in [26l Theorem VI. 2.1.] 
applies only to closed sectorial forms and will thus not be applicable in the present 
case. Instead, we shall borrow results from which generalize the well known 
pseudo-Friedrichs extension [26l Theorem VI. 3. 11.] to quadratic form perturba- 
tions. 

Hypothesis 2.1. Let Hq be a self-adjoint operator in a Hilbert space %, and let 
be a sesquilinear form on T-L such that for all u, v G 2?(t>) 

(2.1) \x>[u,v]\ < \\Hl^u\\\\Hl^v\\, U ,uG27(t>), H a , b := a + b\H \, 

for some a, b > 0, where T>(y>) is a core for llfol 1 ^ 2 - 

Clearly, any sesquilinear form D satisfying the assumptions of Hypothesis 12.11 
may be extended to a sesquilinear form on Q := V^H^ 1 / 2 ) for which l|2.ip con- 
tinues to hold. We will therefore always assume that T>(b) = Q. 

Remark 2.2. Condition (|2.ip is equivalent to 

2\t>[u,v]\ < (a\\u\\ 2 + b\\H \ 1/2 u\\ 2 ) + ^a\\v\\ 2 + b\\H \^ 2 v\\ 2 ) 

If is a symmetric form, then (|2.1|) is equivalent to 

(2.2) Mu]\<a\\u\\ 2 + b\\H \^ 2 u\\ 2 : u G X>(t>). 
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In general, (|2.2[) implies (|2.1[) with a, b replaced by 2a, 2b. 

By the Riesz representation theorem, the formula 

(2.3) (C a , b u,v) = ti[H-l /2 u,H-l /2 v], u,veH, 

defines an operator C a ,b & L(H) with |jC a ,6|| < 1- Assume that the operator-valued 
function C a ,b ■ p{Ho) —> C, 

C a , b (z) := (H a - z)H~l + C a , b , z G p(H a ), 

has nonempty resolvent set, i.e. there exists zq G p(Hq) such that C a; f,(zo) has a 
bounded inverse. By [43l Theorem 2.4] there exists a unique operator H associated 
to the quadratic form f) = f)o + 0, where 

ff [u,v] := (|H o r /2 w,sign(£f )|i/or /2 f) , u,v G Q 

is the form represented by Hq. More precisely, H is given by the formulas 

(2-4) H-z = H 1 a / b 2 C a Az)H 1 J b 2 , 

(2-5) H*-z = H 1 J 2 C a , b (zrH 1 a /2 , zeC, 

and T>(H) is a core for |f/o| • Whenever C a ,b(z) is boundedly invertible, then 
z G z G p(i?*), and 

(2-6) (H - z)- 1 = H-]J 2 C a Az)- l B- a ]J 2 , 

(2-7) (if* - z)- 1 = H- a l /2 d a . b (z)-*H- a l /2 . 

The construction does not depend on a, 6. 

Remark 2.3. If x> is the form of an operator V, 

t}[u, v] = (Vu, v), B,«6 T>(V), 

then H is the pseudo-Friedrichs extension of Hq + V. 

We note that the construction of H in [43] is accompanied by a spectral in- 
clusion, see [131 Theorems 2.11., 2.16., 3.1.]; compare also [TT] for related results. 
For example, [331 Theorems 3.1.] states that if Ho has a spectral gap (—5,S) and 
a + b 5 < 5, then 

(-6 + a + b 5, 6 - a - b 6) + iR C p(H) . 

2.3. Spectral projections. Let S(X — > X) be an operator in a Banach space 
X, and let Q± G L(X) be a pair of complementary projections, i.e. <5+ + Q- = I- 
Then S is said to be decomposed according to 

(2.8) X = Q+X + Q_X 
(compare [5S1 III. 5. 6.]) if 

(2.9) Q±D(S) C V(S), SQ±x = Q±Sx, x G V(S). 
With respect to the decomposition (|2.8j) . 5 is then block-diagonal, 
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where S± := S\q ± x denote the parts of S in Q±X. Clearly, c(5) = cr(S+)Ua(S-), 
and S± are closed, densely defined etc. if and only if S is. We are interested in 
the case where the union is disjoint; in particular, when iR C p{S) and 

(2.10) ff(S ± ) = ff(S)nC ±] C± := {zeC:±Rez>0}. 

If at least one of the the sets a(S) fl C± is bounded, then Q± may be defined by 
the Ricsz-Dunford functional calculus. If S is a sclf-adjoint operator in a Hilbcrt 
space, then this may be accomplished by the self-adjoint functional calculus, even 
if both sets are unbounded. In either case, 

(2.11) — / {S-z)- 1 dzx = Q+x-Q„x, xeX. 

m JxR 

In the general case, the problem of separating the spectrum at infinity arises. The 
following theorem was proved in [29] and is based on JT6l Theorem XV. 3.1.]. We 
slightly simplify the assumptions stated in [29] Theorem 1.1.]; they are equivalent 
by a straightforward Neumann series argument. 

Theorem 2.4. Let S be a closed, densely defined operator in a Banach space X 
such that iR c p(S), limui^oo IKS' — ir/)^ 1 1 = and 

- [ (S-z)- 1 dz 

7" JiR 

exists in the strong operator topology. Then there exist complementary projections 
Q± in X such that S is decomposed according to X = Q+X + Q-X and such that 
(|2T0 ll -(pTT j) hold. 

Remark 2.5. We call Q± the spectral projections corresponding to the right and 
left half planes C±. We remark that S is bisectorial under the stated conditions 
and that the spectral projections may in principle also be defined by the functional 
calculus for such operators, see e.g. [30]. However, the two notions need not 
coincide; in particular, the spectral projections defined by the functional calculus 
may be unbounded. 

We will need the following perturbation result: 

Theorem 2.6. Assume Hvvothesis 12.11 and that (|2.ip holds with a, b > 0, b < 1. 
Then there exists a unique closed densely defined operator H which is associated 
to the quadratic form \) = f)o + 0; moreover, T>(H) is a core for IHq] 1 ^ 2 . If 
iR C p(Ho) n p(H), then the assertions of Theorem \2A\ hold for H . 

Proof. The first part was proved in |43[ Theorem 2.11]. 

To show that H satisfies the assumptions of of Theorem l2.41 we prove that the 
integral 

/oo' 
{(H-b^-iHo-iv)- 1 ) d V> 

exists in the norm operator topology and that 

(2.12) sup|?;| II (H -ir])- 1 - (Ha-iri)- 1 ^! < oo. 

??GE. 
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Since the assumptions of Theorem 12 .41 are obviously satisfied for the self-adjoint 
operator Ho, it then follows that the same holds true for H, by virtue of the identity 

(H - iri)- 1 = (H - ivV 1 + ((H - ir?)" 1 - (H - ty)" 1 ) . 

By the spectral theorem for self-adjoint operators 

(2.13) || (Ho - ir,)- l H a>b \\ < sup 4==> V G R. 

|t|>« v f + V 

A straightforward computation yields for the supremum above (see e.g. [43] ) 

b\t\ 



(2.14) sup "-L^L = rr Va a + ^ 2 . ^R- 
|t|>*v* +if ™ 

Since b < 1, it follows from (|2.13j) and ()2.14j) that for each b G (b, 1) there exists 
R > such that 

(2.15) IKffo-i,,)- 1 ^!! <6, r) 6 R, |r?| > i? 

We may assume without loss of generality that (|2.15|) holds for all rj G R. Other- 
wise, the existence in norm of the two integrals 

((ff-iTj)- 1 -^-^)- 1 ) dr?, / ((tf-i^-Offo-i??)- 1 ) d»j. 

is shown separately. But since the integrand is continuous as function of rj G R, 
the first integral above always exists, while for the second, (|2.15|) holds. 
By (|2.6[) . we have for all rj G II 



(H - ir?)- 1 = Hj /2 ((Ho - ir})H-} + C a , h ) " H Q 



i,6 



= (Ho - irrXi 2 (1 + C a , b (H - ir?)" 1 ^) 1 H^ 2 

oo 

= (Ho - W)" 1 ^ 2 £ [-Ca,6(^o - ir?)" 1 ^,,]" H- 1 / 2 
™=° 

= (H - iT?)- 1 + (^0 - ir?)" 1 ^ 2 E h^, fc (H - i'?)^^]" ^J 7 ' 

n=l 

The sum above converges absolutely by (|2. 15|) and because ||C a ,6|| < 1. Hence, 

oo 

(H - iry)" 1 - (Ho Ir})- 1 = (Ho - ir?)" 1 ^/ 2 E [-C a , fc (H - i^H^f H~^ 2 

n=l 

OO 

(2-16) = (Ho - ir?)" 1 ^/ 2 E ^HCHo - i^" 1 ^/ 2 , 

where 

D n («) := [-C aib (H - z)- 1 ^]" (H - ^H"^ z G p(H ). 
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Note that D n {z) contains n factors of C a ,b and n — 1 factors of (Ho — z) 1 7J a ,b- 
Therefore, by (pT5| . 

\\D n {z)\\ < \\C a , b \\ n \\{H Q - z)- 1 ^^- 1 <b n -\ 

We set 

G a . b (z) := {Ho - z)- x H X Jl; e L(U), z G p(ff ). 
Then, for u,v £ H, 



(([(if-ij?)- 1 -^-^)- 1 ]^^)! 



^ D n (b])G a . b (ir])u, G a , b (iri)i 



(2-17) < (1 - &)- 1 ||Gf a)6 (ir 7 )«|| ||G a , 6 (ii7)t;||. 

Let (/3 ra )neK C (0,oo) be such that p n — > oo, and define {T„} ng ]N by 

T n := 1 f" {(H - i,,)- 1 - (Ho - h?)- 1 ) d V . 
711 J- Pn 

Note that, since the integrand is norm-continuous, the integral exists in norm. 
By (|2.17j) and the Cauchy-Schwarz inequality, we have for all u,v € H, 

" {([(H-ir])- 1 -(Ho-b^^v)] d V 

Ptl 

< (1 -b)- 1 r ||G a , 6 (i77)u|| \\G a ^)v\\ dr, 

J-Pn 

<{l-b)- l ^J ||G a , 6 (ir;H| 2 dr^J (J ||G a , h (ir;) W || 2 dr;J 

= 7r(l-b)- 1 (a\\H - 1 \\+b) \\u\\\\vl 

where the last equality is a consequence of the spectral theorem. Therefore, 
{(T n u, i>)}n£]pj converges uniformly for u,v in the unit ball of T-L. By [551 P- 150], 
it follows that {T n } n€ K converges in norm. 

Another application of the spectral theorem yields, using (|2 . 1 T[) . 

UH-b!)- 1 -(Ho-iv)- 1 ^ sup \([(H -irj)- 1 -(Ho-ivr^^v^ 

IM|=|MI=i 

^:„M|2 /-, 7,-1 „..„ a + b\t\ 1 



< (l-br'WGa^)^ = (l-b)- L SUP 



t& r&o) H+tM~ l \n\ 
<(i-br - P i+M i <(i_6)-^(f||^|| + 6) * 

te<r(i? ) 2 I*I I 7 ?! V2 7 I 7 ?! 

This proves (|2T2|) . □ 
2.4. Graph subspaces and angular operators. 

Definition 2.7. Let X, Y be Banach spaces and Z = X + Y. A subspace C C Z is 
called a grap/i subspace with respect to X if there exists an operator A x € L(X, Y) 
such that 

£ = {x + A x x : x e X} . 
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In this case, Ax is called the angular operator of C with respect to X. 

For simplicity, we shall also call a subspace M. C Y of the form 

M ={A Y y + y:yeY}, 

Ay G L(F, X), a graph subspace, although the term inverse graph subspace" would 
be more appropriate. 

Remark 2.8. Let Px be the projection of Z onto X along Y and Py = Iz — X. 
It is easy to see that C C Z is a graph subspace with respect to X if and only if 

Px\c-C^X 

is an isomorphism and that the angular operator Ax is given by 

Ax = Py{Px\c)- 1 - 

Proposition 2.9. Let X, Y be Banach spaces, Z = X + Y , and let £, M. C Z be 

graph subspaces with respect to X and Y, with angular operators Ax € L(X,Y) 
and Ay £ L(Y, X), respectively. Then the following are equivalent. 

i) Ix — Ay Ax has a bounded inverse; 

ii) Iy — Ax Ay has a bounded inverse; 

iii) The operator 

w -=(Z t) EL(z) 

has a bounded inverse, with 

iv) C + M = Z. 

Proof. The equivalence of i)-ih) follows from the Schur-Frobenius factorization; for 
example 

'IX \ fix Ay 

v Iy - AX Ay J \ Iy 

The formula (|2.18[) is easily verified by a direct computation. To prove the equiv- 
alence of iii) and iv) , we observe that since 

'x\ = ( x \ + ( Ayy^ 
J/J \A x x) \ y 

W is surjective if and only if C + M = H, and W is injective if and only if 
CD A4 =0. An application of the closed graph theorem completes the proof. □ 

In the Hilbert space situation, we have the following useful proposition see e.g. 
(23 Corollary 3.4]. 

Proposition 2.10. Let H be a Hilbert space and C, M. C % subspaces with cor- 
responding orthogonal projections P^, Pm- Then \\Pc ~ Pm\\ < 1 if and only if 
Ai is the graph of an operator K £ L(£, C ). In this case 

(2-19) \\K\\= J P f- PM } unt \\P £ -P M \> im 



Ix 


Ay\ 


fix 




\A X 


Iy 


\Ax 





w 



^1-\\Pc-Pm\\ 2 V^+W\ 
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It is easy to see that d(C,M) := \\Pc — Pm\\ defines a metric on the set of 
subspaces of a Hilbert space H. Let us introduce the angular metric by 

d a (C,M) := arcsin||Pc - Pm||- 

The fact that d a is indeed a metric (i.e. satisfies the triangle inequality) was proven 
in [5]. The angular metric is related to the operator angle 



@c,m : = arcsin \/lc~ PcPm 
between the subspaces C and M. by the equality d a {C,M) = ||0£,jw||- 

2.5. Accretive operators in indefinite inner product spaces. Our main 
tool in the proof of Theorems 13.11 and 13.31 is the following result about accretive 
operators in indefinite inner product spaces. It is a slight generalization of |29[ 
Theorem 1.4.]. 

Definition 2.11. Let TL be a Hilbert space and W a bounded self-adjoint operator 
on H. 

i) An operator T(TL — > H) is called W-accretive if 

Re {WTx, x) > xeV(T). 

ii) A linear manifold C C % is called W-nonnegative (W-nonpositive) if 

(Wx,x) > ( < 0), x e C. 

Theorem 2.12. Let % be a Hilbert space and W a bounded self-adjoint operator 
on %. Let T(H — > %) be a closed, W-accretive operator such that iR\ {0} C p(T). 
Assume that the integral 

(2.20) ^ / (T-z)~ l dz, 

m JiR 

exists in the weak operator topology and is the difference of two complementary 
projections Q± G L('H), 

(2.21) —[ (T-z)- 1 dz = Q+-Q_. 

7" JiR 

Then Q+H C % is a W-nonnegative subspace and QJH, C % is a W-nonpositive 
subspace. 

If £ p{W), we denote by P± the spectral projections onto the positive and 
negative spectral subspace ofW. Set H± := P±H and W± := ±Vl^|^ ± . Then there 
exist operators K± € L('H±,'H = p) such that 

(2.22) Q+H = {x + K + x:xeP + H}, QJH = {y + K_y : y g P_H} . 
Moreover, 
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Proof. Let x e Q+H. Then 

[x,x] = Re[x, x] = Rc[(Q + — QJ)x,x] 



Re f — / [(T - z)~ 1 x,x}dz ) 
\ m Jin J 

-[ Re[T{T-iT 1 )- 1 x,{T-r ] )- 1 x}dT]>0. 
n Jn 



Thus, Q+H is nonnegative. Analogously, one shows that Q-H is nonpositive. 

If € p(W), then H equipped with the indefinite inner product (W-, •) is a Krein 
space with fundamental decomposition H = H+[+]H-. Since Q+ + Q- = I, 
[TJ 1.1.25] implies that the subspace Q+H is maximal nonnegative and Q+'H is 
maximal nonpositive. It follows that there exist operators K± € L('H±,'H^) 
such that ([2~2^|) holds, see e.g. [1 Theorem 11.11. 7] or [2S|. Moreover, K± are 

1 /2 

contractions with respect to the canonical norms | W± ■ \\ in P±"H induced by VK. 
Stated differently, 

\\wt 1/2 k ± w±~ 1/2 \\ < 1, 

which implies 

\\K ± \\<\\w-^nwn^-^^ 



y inf <j(W t ) ' 

The last equality is a consequence of the spectral theorem for self-adjoint operators. 

□ 

More advanced results of this kind, stated in terms of interpolation spaces, may 
be found in [31] . 

3. Main results 

Theorem 3.1. Let Hq be a self-adjoint operator with a spectral gap (—6,6) in a 
Hilbert space 7i. Let be a quadratic form such that T>(t) is a core for |iJ | 1//2 
and such that 



\H \-^xAH \-^ 



< 1, 



sup 

IMI=lly|l=i 

H Q \-^x,\H Q \-V*J 



sup 

x.yGP+HUP.H 

IMHMI=i 



=:P< 2- 



Furthermore, let P± be a pair of complementary projections in % such that 



v := \\P± — P±\\ < 1, arctan . / — h arcsin^ < — . 

y 2 — 3p 2 

Then the following hold: 

i) There exists a unique closed densely defined operator H associated to the 
quadratic form f) := f)o + 0, and T>(H) is a core for IHq] 1 ^ 2 . 
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ii) There exist complementary projections Q± in H such that H is decomposed 
according to H = Q+H + Q-H , and 

<j(H\ Q+n )=a(H)n€ ± . 

hi) The restrictions P±\q ± u '■ Q±H ~^ P±H are isomorphisms, and, with 
X± := P T (P±\Q ± ny 1 , we have 

Q+H = + X+x : x e P+h} , Q-H = \y + X.y : y e P-H^ , 
||X ± || < tan ( e 



arctan 



2-3p 



arcsm v 



Remark 3.2. It is easy to see that for any invariant set M C H of Ho, the 
following are equivalent: 

i) There exists u> > such that 



sup 

x,y£M 
\\x\\ = \\v\\ = l 



|tfo|- 1/2 *,|ffo|- 1/2 y 



< w; 



ii) There exist a, & > 0, a + & £ < such that (|2 . 1 1) holds on M. 

Theorem 3.3. Let H$ be a self-adjoint operator with a spectral gap (—(5, S) in a 
Hilbert space H. Let be a symmetric quadratic form such that 2?(t>) is a core for 
jiJol 1 ^ 2 and such that 



sup 

11*11 = 112/11 = 1 

sup 

IMI=IMI=i 



\H \-^x,\H \-^y 



\H \-^x,\H \-^y 



< 00, 



P<1. 



Furthermore, let P± be a pair of complementary projections in H such that 

n 



v 



\P± — P±\\ < 1, arctan 



p 

h arcsm v < — . 

2- p 2 



Then the following hold: 

i) There exists a unique self-adjoint operator H associated to the quadratic 
form f) := f)o + moreover, € p(H), and T){H) is a core for \Hq\ x > 2 . 

ii) Let Q± denote the spectral projections of H corresponding to the positive 
and negative spectrum, respectively. Then P±\q ± u : Q±H — > P±H are 
isomorphisms, and, with X± := P+{P±\q ± u)~ 1 , we 



)+H = {re + X+x : x e P+h} 
\\X±\\ < tan ( arctan , 



2-p 



M = {y + X_y:yeP-U) 
arcsm// I , X- = —XI. 
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Remark 3.4. If P± are orthogonal projections, then the result of Theorem 
may be equivalently stated as 

\\P± — Q±\\ < sin ^arcsin + arcsin v 

An immediate consequence of Theorems 13.11 and 13.31 is a block-diagonalization 
of If. 

Corollary 3.5. With respect to the decomposition % = P+% + P—H, we have 

i x.V 1 „ ( i x\ ( X o 



(3 - 1} [x + I ) H yx + 

where Z± are similar to H\q ± u. 

If v = and (H — i?;) -1 — (Hq — it])^ 1 G iSp('H) for some (and hence for al 
T) S R, then X± G S P (H). 

Proof. We set 

Since Q±H are complementary graph subspaces with angular operators X±, it 
is easily seen that Q±W = WP± and that W is bijective and hence boundcdly 
invertible by the closed graph theorem. Since H is decomposed according to 
Q+T-L + Q-H, it follows that W~ 1 HW is decomposed according to P + H + P-H, 
i.e. it is block-diagonal with respect to this decomposition, and (|3.1j) holds for 
some closed densely defined operators Z± on P±H. Noting that 

W\p ±n = (P ± \ Q±n r 1 , 

we find that 

Z± = {P±\ Q± u)H\ Q± u{P ± \ Q± u)- 1 . 
If the resolvent difference of Hq and H belongs to S p (7i), then 

Q±-P±= n-lim |i f (H- ir))- 1 - (H - irj)- 1 d v e S p {%). 
p^oo 2tt J_ p 

With respect to the decomposition % = P+W + P-H, e.g. the projections P + 
and Q + are given by 

r*>>\ p -f 1 °\ n -( V-X-X+r 1 -(i-x-x+y^- 
[6 - Z) ^+~\o oj' y+ ~ \X+(I - X-X+)- 1 -X+(I -x_x+)- l X- 

see [U p. 63]. It follows that 

x+ii-x-x+y 1 es p (u), 

and hence 

x+ = x + (i - X-X+y 1 ^ - X-X+) es p {u) □. 
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Remark 3.6. In the case when t> is symmetric, the fact that X- = XI implies 
that the bounded operators 1 — X^X± are self-adjoint and uniformly positive. The 
transformation W in Corollary [33] can then easily be made unitary by multiplying 
from the right with diag(fi+, fi_), 

n± := (i-x T x ± y 1/2 , 

see [55]. In fact, the resulting operator takes the more familiar form (on the 
right-hand side) 

(3.3) X Q n -^j =[I- (Q+ - P + f}~^[P + Q + + P.Q-], 

as can be checked by a straightforward computation, using Phroposition [^TTUl and 
(|3.2p . This is the direct rotation between the subspaces P±H and Q±H, see [15] . 

For the next two theorems, we restrict ourselves to the case P± = P± . The aim 
is to approximate the exact block-diagonal operator by simpler operators which 
are amendable to computation. To this end we introduce a coupling constant 7 
for the perturbation and consider the family of operators 

H(rf) := H + 70, 7 GB, 

where D C C is the open unit disk. Notice that Theorem 13.11 i) is valid for the 
whole family H(-f), 7 G ID, while ii) only holds for 7 G (—1,1) since 7D is not 
symmetric for 7 G C\l. Therefore, under the assumptions of Theorem 13. H i), the 
following operators arc well-defined for any 7 G D: 

• The spectral projections Q±(j) of #(7) corresponding to C±; 

• The angular operators X±( 7 ) := P T (-P±|q ± ( 7 )) ; 

• The inverse square roots fl±( r y) := (1 — X T ( , -f)X±('-f)) ' . 

The latter can no longer be defined by the spectral theorem. Rather, since 
sup 76D ||^"±(7)|| < 1, we may define them by a norm-convergent power series 

OO y . \ 

(3.4) n±( 7 ) = V ( ~T ) (-x T ( 7 )x ± ( 7 )f . 



Let Hdiag(7) be the block-diagonal operator (with respect to the decomposition 

n = p+n © P-U) 



iW7) : = ^(7)^ (7)^(7), W( 7 ) 
Moreover, if is symmetric, define i?diag(7) by 
£dia g (7) := C/( 7 )- 1 H(7)f/(7), Uin) : 



/ X_( 7 ) 



n+( 7 ) X_( 7 )Q_( 7 ) 

x + ( 1 )n + ( 1 ) n_( 7 ) 

We show that the operators H^ iag (j), (j) corresponding to the formal iV-th 

order Taylor polynomials of i?diag(7), respectively -ffdiag(7), converge in the norm- 
resolvent sense to the exact block-diagonal operators as N — > 00. Furthermore, we 
make precise in what sense these operators may be viewed as Taylor polynomials. 
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For technical reasons we assume that is the form of an operator V, i.e. 
v[u,v] = (Vu,v), u,v G X»(t)) = T>(V), 

Theorem 3.7. Assume that t) in Theorem 13.11 is the form of an HQ-bounded 
operator V with H^-bound less than [2 — 3p)/(2 — 2p), then {-ffdiag(7)}76D is 
a holomorphic family of type (A) with 2?(-ffdiag(7)) = ?)(Ho). For 7 € ID and N 
large enough, the operators H^ iag (j), defined by 

N -y" d" 

^dLg(7)« ~ E ^1 d^T (^di-g(7)«) (7=0, u € 2?(ffo), 

n=0 ' ' 

Ziawe nonempty resolvent set and H^ iag (-f) — > -ffdiag(7) as N ^ 00 in the norm- 
resolvent topology, uniformly on compact subsets ofD. 



Theorem 3.8. Assume that t> in Theorem 13.11 is the form of an HQ-bounded 
symmetric operator V such that \\VH ~ 1 \\ < 1. Then {ifdiag(7)}7eB is a self- 
adjoint holomorphic family, and 



T( 7 ) := |H |-i/2tf diag ( 7 )|tf |-i/2 ; 7 G D, 

is bounded-holomorphic. For 7 G (— 1, 1) and N large enough, there exists a unique 
self-adjoint operator H^^j) associated to the quadratic form 

N <v n H™ / \ 

ftSL«(7)[«,«] -E^r^r (TtTjiffoi^Viflbi^t,) 



n=0 



7=0 



and X>(^ ag (7)) zs a core /or liJol^ 2 - Moreover, H^ mg (j) -> #dia g (7) as TV ^ 00 
in i/ie norm-resolvent topology, uniformly on compact subsets of (— 1, 1). 

Remark 3.9. The norm-resolvent convergence in Theorem 13.71 implies conver- 
gence of the spectra and spectral projections of H^^) to those of #diag(7) since 
all operators are self-adjoint, see [3H VIII. 7.]. 

On the other hand, in Theorem 13.11 norm-resolvent convergence only implies 
that if z € /o(-ffdia g ( 7 )), then z e p{H^ iag {j)) for N sufficiently large. By con- 
traposition, this means that if z G <T (^di a g(7)) f° r a ^ N sufficiently large, then 
z e cr(iJdia g (7))- However, it can easily be seen from the proofs of Theorems 12.61 
and 13.71 that the spectral projections corresponding to C± do converge. 

Remark 3.10. Like in [37], the method to prove convergence is based on analyt- 
icity. We remark that even if is symmetric, the detour through non self-adjoint 
operators (i.e. complex 7) is unavoidable. Indeed, in order to estimate the radius 
of convergence of the polynomial approximation, we need complex analyticity, not 
just real analyticity. 

4. Application to Dirac operators 

4.1. The free Dirac operator. Let us denote the free (i.e <f> = 0, A = 0) Dirac 
operator defined on the first order Sobolev space by Hq, i.e. 

(4.1) ff =(_ i j. v _1 °' 1 V )' V(H ) = H 1 (R 3 ,G 2 )®H 1 (R 3 X 2 ). 
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It is well-known that that the spectrum of Hq is absolutely continuous and equal 
to 

(4.2) cr(# ) = (-00, -l]u[l,oo). 

The Foldy-Wouthuysen transformation J7fWi m momentum space a multiplica- 
tion operator by the function 

(4.3) u( P ) := L±m±frL* -^_(i+m -p 



N(p N{p) \ -a ■ p 1 + E{p) 

where 



E(p) := sjl+p\ N(p) := y/2E(p)(l + E(p)), p := |p|, p 6 R 3 , 
block-diagonalizes the free Dirac operator: 



/ -ier-V\ /VT^A 

UFW 1 :_ v7 r I ^FW 



-icr-v -/ y * w V -VT^A 

Moreover, the orthogonal projections P± = 1/2(1 ± i? |ffo| _1 ) onto the positive 
and negative spectral subspaces of Ho in momentum space are given by 

{ ' A±[P) ~2E(p){ <r-p Tl + E{p) 

For fixed p € R 3 , the matrices A±(p) are orthogonal projections in C 4 , and 
both possess a two-dimensional eigenspace corresponding to the eigenvalue 1. For 
example, given two orthonormal vectors ui, it 2 € C 2 , the vectors 

yb-- -J- f( 1 + ^M i-1 2 

form an orthonormal basis of the eigenspace belonging to A + (p). It is thus seen 
that 

P,U - ( (1 - + £( i! S - H ) 



«(.) v <r- (■)»(•) 7 •"="»/ ~ \V- ¥ +" 

where T is the Fourier transformation, u := J 7 !* and 
-1 *■(■ 



A+ : = j- 1 ^— -Ai-JF= -i.<7. v (i + a/i - Aj eL(H u ,H0 
is the angular operator of P+H with respect to ?i u . 

4.2. Assumptions on the potential. We impose the following general condi- 
tions on the potential V. 

Hypothesis 4.1. Let V be a measurable Ax 4-matrix-valued function, and assume 
that there exist a, b > such that 

b 

(4.5) II^MIluc 4 ) < a + 1 — fi for almost all x 6E . 
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In most physically relevant applications, V is given by 

However, we will not assume V to be of that particular form. Most importantly, 
our assumptions cover the Coulomb potential V = Za/\ ■ |, where a ~ 1/137 is 
the fine structure constant. 

The following proposition summarizes some useful inequalities. 

Proposition 4.2. Let Hq be the free Dirac operator (|4.1j) . Then the following 
hold. 

i) II I • < 2 \\H u\\ for all u 6 V{H Q ); 

ii) |(| • r 1 ^,^)! < tt/2 II |#o| 1/2 «|| 2 /or we Q; 

hi) |(| • |-Vu)| < 1/2 • (tt/2 + 2/tt) I) l^ol 172 "!! 2 V nHtieP+fiUP.fi. 

Inequalities i) and ii) are the classical Hardy, respectively Kato inequalities. 
Inequality iii) is due to (40l [41] and [6] , see also [2] . It establishes the boundedness 
from below of the Brown-Ravcnhall operator and was also used in |13j to prove 
variational principles for Dirac operators with Coulomb-like potentials. 

4.3. Magnetic fields. If we want to include magnetic fields (i.e. A ^ in (|4.6j) ). 
then the condition (|4.5|) is too restrictive since it does not allow e.g. for constant 
fields B = curl A. The magnetic field term should thus not be considered as a 
perturbation; rather, we should include it in the definition of the unperturbed 
operator. If, for example, A e L 3 oc (]R 3 , K 3 ), then the magnetic Dirac operator 

is essentially self-adjoint on Cq°(M? \ {0}) 4 , see e.g. [3£l notes sect. 4.3]. In the 
following, we denote its self- adjoint closure by the same symbol Owing to 

the special structure of the operator matrix, (—1,1) is a spectral gap for Hj^. 
By combining Hardy's and Kato's inequalities with the diamagnctic inequality for 
Schrodinger operators [5J Theorem 4.5.1.], one can prove the following proposition, 
see e.g. [23l (4.7), (4.9)]. 

Proposition 4.3. Assume that A <G i 3 oc (]R 3 , K 3 ) and that B := curl A is essen- 
tially bounded. Then the following hold. 

i) II I • r^f < 4 ||ff A w|| 2 +4 • llBHoo |M| 2 for all u G V(H A ); 

ii) KH-Vu)! < 7 r/2|||F A | 1 /\|| 2 + 7 r/2.yp|Ulk||Vora^neX'(|HA| 1/2 ). 

As was noted in [24] , the boundedness assumption on the B-field may be relaxed 
by using an estimate due to Balinsky, Evans and Lewis [3] , relating the Schrodinger 
operator (V — iA) 2 to the Pauli operator (or • Va) 2 , when the latter has no zero 
modes. 
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Proposition 4.4. Assume that A G L 3 oc (R 3 , R 3 ), and let B := curl A 



-1/2 



(4.8) <5(B):= inf ||(1 - S*S)x\\, S := |B|^ ( ((7 . y A ) 2 + |B|) 

xGH, I a; II = 1 

T/ien < (5(B) < 1, and the following hold. 

i) || | • < 2 ■ 5(B)" 1 \\H A uf for all u G Z>(ff A ); 

ii) |(| ■ I -1 !*,?/)! < vr/2 ■ (5(B)" 1 ||i? A w|| 2 for all u G P^aI 1 / 2 ). 

4.4. Block-diagonalization and convergence of the DKH method. For no- 

tational convenience, we subsequently identify the abstract unperturbed operator 
Hq in Section 3 either with the free Dirac operator (denoted by the same symbol 
above) or with the magnetic Dirac operator Hj±, depending on whether A vanishes 
or not. Correspondingly, the projections P± in Section 3 are identified with the 
spectral projections of the free or the magnetic Dirac operator. 

Proposition 4.5. We have \\P U - P+\\ = \\Pi - P-\\ = 1/V2- 

Proof. For the free Dirac operator, this can easily be inferred from (|4.4|) . For the 
magnetic Dirac operator, it follows from supcrsymmetry arguments, see [29]. □ 

In view of Remark 13.21 and Propositions 14.2114.41 sufficient conditions for The- 
orems 13.11 13.31 13.71 13.81 to hold may easily be obtained in terms of a and b. For 
example, by Proposition 14. 5[ in the absence of magnetic fields (A = 0), the as- 
sumptions of Theorem 13.31 arc satisfied for P + = P u , P_ = Pi (i.e. for the natural 
decomposition (|1.2|l ) whenever 

b (tx 2' 
2 + 



<1, 



whereas for Theorem 13. 81 the left hand side has to be less than 1/2. In particular, 
for the Coulomb potential V — Za/\ ■ \ the above inequalities amount to Z < 124 
and Z < 62, respectively. In the case A ^ 0, the upper bound for Z depends 
on the magnetic field; Proposition 14.41 yields Z ■ (5(B)" 1 < 87 (as compared to 
Z-SiB)- 1 < 60 in [24]). 

5. Proofs of the main results 

Proof of Theorem 13.11 i) follows from [331 Theorem 3.1] and Remark 
addition, we have iR C p{H). 

ii) is then a direct consequence of Theorem 12.61 



iii) is shown in two steps. First for v = and then for the general case. 
We claim that H is WVaccretive for i = 1, 2 and 

^ : =2^ P +- P - W >--= P +-2hp P - 

Theorem l2 . 1 21 then implies the existence of the angular operators X± of Q±H with 
respect to P±H such that 



2-3p 
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This completes the proof in the case v = 0. 

For the general case, let P± and Q± be the orthogonal projections onto P±H 
and Q±H, respectively. By Theorem 1.6.35.], we have 



Hence, by the triangle inequality for the angular metric and Proposition 12 . 101 
arcsin \\P± — Q±\\ < arcsin \\P± — Q±\\ + arcsin \\P± — P±\\ 
< arcsin ||P± — Q±\\ + arcsin \\P± — P±\\ 



By Proposition 12.101 again, Q±H = Q±H is a graph subspace with respect to 
P±H = P±H, with angular operators X± satisfying 



We now show that H is W^-accretive for i = 1, 2. To this end, for fi± > 0, we 
set W := M+-P+ — H—P-. Since P±Q C Q, the following sesquilinear forms are 
well-defined on Q: 

6d[m, v] := fi+ v[P+u, P+v] — jtx_ o[P_u, P-v], 

6q[w, u] := t)[P-ii, P+u] - /U- t)[P+u, P-u], 



X±[«,«] :=M±ll^o| 1/2 P±«!l !|Po| 1/2 P ± i'|| +n±BcX> u [P±u,P±v]. 
Then, for u g T>(H), 

Re (WfTu, u) = Re (.Hit, Wu) = Refj[u, Wu] = Re f) [u, Wu] + Re x>[u, Wu] 
^ 5 ' 1 ^ = RefjoK W^"] + Rc6dM + Rc6oM = X+N + X-M + Re&oM- 
By assumption, we have 

|Rc6dM| <p-n+\\ Hl /2 P +U \\ 2 +p-fi_ || P 1/2 P-«|| 2 , 

|Re6oM| < <t\Rcv{P-u,P + u}\ <p-^\\ Hl /2 P + uf +p ■ °~ || E^P-uf, 



||P± - P±|| < ||P± - P±||, ||Q± - P±|| < ||Q± - P±||- 




which is equivalent to 





where 



a 



It follows that 



X±M>M ± (i-p)lll#o| 1/2 P ± w|| 2 >0, 
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and, by (|S~T|) . 

M"^") * (* - 2, + 7l- P ) ) ^ + i 1 - 2 /'(I -p) ) 

For //+ = p/{2 — 2>p) and /x_ = 1, the expressions in brackets are nonnegative. The 
same is true for p + = 1 and p_ = p/{2 — 3p). □ 

Proof of Theorem i) follows from gH Theorem 2.16] and Remark [3^1 

ii) The proof is identical to the one of Theorem l3.1l iii) above; the only difference 
is that here we can choose 

a := \p+ - |U_|, 

and H will be W^-accretive for 

W 1 :=-^—P + -P-, W 2 :=P+--^P-. 
2 — p 2 — p 

Since H is self-adjoint, the equality A_ = XI follows from the orthogonality of 
Q+H and Q-H. □ 

Before we proceed to the proof of Theorem l3.7[ let us briefly recall the following 
proposition, which we shall subsequently use without further mention. The proof 
is a straightforward application of the Neumann series, see e.g. [UJ Lemma 1.8.1]. 

Proposition 5.1. Let X be a Banach space and Q C C open. Assume that the 
operator-valued function T(-) : — > X is holomorphic in z$ £ and that T(zq) 
has a bounded inverse. Then T(z) has a bounded inverse in a neighbourhood of z§, 
and T(-) _1 is holomorphic in Zq. 

Proof of Theorem 13.71 We first show that W(-f) is holomorphic in D. 
Replacing C by jC in (|2.16|) . one observes that 

(F( 7 ) - ir?) 1 - (Ho - iv)' 1 

is a norm-convergent power series for 7 G D. By the proof of Lemma 12.61 the 
integral 

Q±{l)-P± = ^ I ((H^-z^-iHo-z)- 1 ) dz 

converges in the norm operator topology, uniformly in 7 on compact subsets of D. 
This shows that {i7(7)} 7g B and {Q±(7)} 7 eB are holomorphic families. 
From (13.21) we infer that 



*+(7) =g 2 i(7)Qii(7)" 1 , 7€B, Q+h)=-(Qij(7))l j=1 , 

whence {X + (7)} 7S ]d is holomorphic. The proof for {A_(7)} 7e B is analogous. 
Consequently, {VF(7)} 7£ d is holomorphic. 

Since the TJo-bound of V is less than one, the family {i?(7)} 7e D is holomorphic 
of type (A), with V{H{^)) = V{H ) for all 7 e D [56] VII]. We now show that 
i?diag(7) is also holomorphic of type (A). It is sufficient to show that 

(5.2) W^MHo) C V(H ), 7€»- 
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Indeed, since Hq is closed and W( 7 ) is bounded, the operator 

F( 7 ) := HoW^H^ 1 
is closed and thus bounded by the closed graph theorem. It follows that 

(5.3) V(H dias ( 1 ))=V(H W(j))=V(Y(j)Ho) = V(H ), 7 e D. 

It remains to show that -ffdiag(7)u is holomorphic for every u G T>(Hq). Since 
H^Hq 1 is bounded and holomorphic in norm, the same applies to -Hdiag^Po -1 
in virtue of 

In particular, Hdi ag (j)u = Pdiag( 7 )Po~ {Hqu) is holomorphic for every u G T>{Hq). 
We now prove (|5.2[) : note that this is equivalent to 

(5.4) X ± ( 7 )P±X>(#o) C P T V(H ), 7 gD. 

Since P±T>(Hq) C T>(Ho), we can write -ff( 7 ) as an operator matrix with respect 
to the decomposition "H = P+H © P_"H as follows: 

_ / A P\ ._ [P + (H + 1 V)P+ P + (H + >yV)P-\ 
a ™ \C Dj- {P_(H + 7 V)P+ P.(H + 7 V)P-J ■ 

By the Schur-Frobenius factorization (see e.g. [42] for unbounded operators), the 
bounded invertibility of (P( 7 ) — b]) is equivalent to the bounded invertibility of 
e.g. the first Schur complement 

S+(iri) ~A- i v - B(D - ir^C, V(S+(i V )) := P+V(H ). 

Moreover, we have 

'S+iiv)- 1 -S+{iri)- 1 B{D-i n )- 1> 



(5.5) (tf( 7 ) - iri)- 1 = 



where the lower entries of the matrix are bounded operators which can be expressed 
in terms of B, C, D and 5+(ir/) _1 ; we won't need the explicit expressions here. 

Since Q+P( 7 ) C H(j)Q+, it follows that Q + {H (7) — ir?) 1 = (#(7) -h/)- 1 Q+. 
Comparing the ranges of the operators on the left and right hand side yields 

Q+H n V(H Q ) = (#(7) - l^Q+U. 



Projecting onto P+'H on either side, we obtain, using 
(5.6) 

{x G P+H : X + (j)x G P-H} = {S , + (i77)- 1 (l - P(P - ijjJ-^.X'+frJa: : a; G P+H} . 
Setting 

p := \\P+VP-(P-H P- -i77) _1 ||, 

q:= ||P_FP_(P_ J ff P_-i?7)- 1 ||, 

r-UFCHo-ir?)- 1 !!, 

and observing that p 2 + q 2 < r 2 , we obtain, by a Neumann series argument, 

\\B{D-i V )- l \\<-^<^==. 

l-q VI - ?' 2 
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Denoting by bo the -f/o-bound of V, we find that 

Hence, for \rj\ sufficiently large, 1 — B(D — h])^ 1 X + (j) is a isomorphism in P+H, 
and the right hand side of (|5.6j) equals P+T>(Ho). Therefore, we have 

{x G P+U : X+(-y)x G P-H} = P+H, 



which is equivalent to the first inclusion in (|5.4[) . The second inclusion is shown 
analogously, by using the second Schur complement 

S-Q.7]) := D - irj - C(A- irfj^B, V(S-(iri)) := P-V(H a ). 

The fact that -ffdiag( 7 ) is holomorphic of type (A) now implies the norm- 
resolvent convergence of its Taylor series. We notice that the proof given here 
can be adapted to the case when Hq is not boundedly invertible (or even when 
H is only a Banach space) by regarding the following operators as maps from the 
Banach space T>(Hq) (with the graph norm) into T-L. 

For N G M and 7 GB, define the operators R N (j) by 

R N ( 7 )u := ffdiag(7)^ - ffdiag(7K u G 2>(fl"( 7 )) := V(H ). 

Since H^^^Hq 1 is bounded- holomorphic in D, its Taylor series converges uni- 
formly on every compact subset K C D, which means that 

sup ||i? Ar ( 7 )F ( 7 1 || 0, iV^oo. 

By the stability of bounded invertibility [5^1 Theorem IV. 1.16.] it thus follows that 
-^diag(7) nas a bounded inverse for sufficiently large N; moreover, by the second 
resolvent identity, 

||ffdia g (7)- 1 - H^-'W < II^Ct)" 1 !! ||i? JV (7)^dia g (7)- 1 || 

< \\H Aiae (j)-i\\ ||^( 7 )g diag (7)- 1 || 

1 - ||^( T )jy diag ( T )-i|| 

The latter converges to zero as N —> 00 since 

||i^(7)^di«(7) _1 || = ||i? JV (7)^ 1 ^(7)ffo^(7)- 1 W r (7)- 1 || 

< 11^(7)^' II IWrtll WHoH^y'w wwfr)- 1 ]]. 

The convergence is uniform in 7 G K since the functions 1 1 -^diag ( ' ) _ 1 1 U 11^(011 e ^ c - 
are continuous from K to [0,oo) and hence take their maximum on the compact 
set K. □ 

Proof of Theorem 13.81 Like in the proof of Theorem 13 .71 one shows that K± (7) 
are holomorphic in D. The absolute convergence of the series 13^41 locally uniformly 
in 7 G D, implies the holomorphy of ^±(7) and hence of {7(7). The equation 

t diag (7)- 1 = [/( 7 )- 1 ^( 7 )- 1 t/( 7 ), 7 G », 
then shows that -ffdia g (7) is holomorphic in D. 
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Clearly, since V is symmetric, {-ff (7)} 7 gB is a self-adjoint family, i.e. 

H(-y)=H*(j), 7 eD. 
By the identity theorem for bounded-holomorphic functions, we have 

al( 7 ) = -x+w, u{i)u^y = uffiufr) = i 

for all 7 £ D. It follows that -ffdiag( 7 ) is a self-adjoint family. 

We now show that for all 7 £ D, |.ffo| _1//2 -ff(7)|-ffo| _1//2 extends to a bounded 
operator by closure. We claim that it is sufficient to show the following: 

(5.7) I±(#±SC Q, || |tf | 1/2 X ± (7)|ffor 1/2 || < L 

Indeed, it then follows that for all u £ P±H, v £ P±Q, \\u\\ = \\v\\ = 1, 

|(r! ± ( 7 )ii/or i/2 ",i^or /2 «) 



E 

n=0 



-1/2 



-X T ( 7 )X ± (7)]"| J ffo|" 1/2 ^,|^o| 1/2 « 



Ef y 2 )(|^r /2 [-X T ( 7 )X ± (7)]"|^ 



n=0 
00 

^E 

n=0 



I 1/2 u, U 



-1/2 



|i/ | 1/2 X_(7)| J ff r 1 / 2 \H Q \ 1/2 X + (j)\H Q 



-1/2 



By the definition of the adjoint, this implies that 

^±(l)\H r 1/2 P±HcP ± Q, 
whence, by the closed graph theorem, 

\H \ 1/2 n ± ( 7 )\H \-^ 2 £L(P ± H). 
It is then easy to see that 

(5.8) |H | 1 / 2 (7(7)|i/or 1/2 e HH), \H Q \^ 2 U^y\H \-^ 2 £ L(H), 
and we have 

i^or 1/2 J ffd iag (7)iffor 1/2 



H \ 



/ 2 C/(7)1^ | 1/2 ) (|H r 1/2 H( 7 )|i?or 1/2 ) (|i/ | 1 / 2 C/(7)|ffo|- 1/2 



Since the first and the third factor above are bounded, the claim is proved if the 
second factor has a bounded closure. Since T>(H(-f j) = V(H ). it follows by the 
Heinz inequality that X>(|i/(7)| 1/2 ) = X>(|ff | 1/2 )- The latter is equivalent to the 
boundedness of \H \-^ 2 H ( 7 )|# r 1/2 on \H \ 1/2 T>{H{j)), see e.g. H3 Theorem 
3.2.]. Since if (7) is a core for l-f/ol 1 / 2 , this domain in dense in "H, and the operator 
has a bounded closure. This may also be verified directly by formula (|2.4p . 



To prove (|5.7p . we introduce the operator 

W := Li\H \- 1 P + -\H \- 1 P_ £L(H), ,i 



\VH n 



< 1. 
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For u <G T>(Hq), we have 

Rc(WH(-f)u, u) = Rc(iJ( 7 )u, Wu) > (H u, Wu) - \{Vu, Wu)\ 

> /i||P+w|| 2 + \\P-u\\ 2 - \\u\\ \\VH^\\ \\nP + u - P-u\\ 

> a \\P + uf + ||P_ U || 2 - \ ■ \\VH^\\ {\\u\\ 2 + M 2 ||P +U || 2 + ||P_ U || 2 ) 

> (m - \ ■ \\VH^\\ - ix ■ \ ■ \\VH^ \\P + uf + (1 - WVH^W) \\P_u\\\ 

and both summands in the last line are nonnegative. Thus, -ff (7) is VF-accretive, 
and Q±{j)V. are VF-nonnegative and VF-nonpositive, respectively, by Theorem l2.12l 
From the PF-nonncgativity of Q + e.g. it follows that for u € P+H, 

,4 m-^\\> - \w-^x M uf = (w ( x+ » 7)o ) . ( x+ » 7)o )) > o, 

or, put differently, 

|| |tfor 1/2 * + (7)l#o| 1/2 u|| < f,\\u\\, u e P+V(H ). 

Hence, |i/ | _1 ^ 2 ^+(7)|P^o| 1 ^ 2 has an extension to an operator in L(P + 'H) bounded 
by fi, and, by duality, 

ip i i/2 * + (7)*i#or i/2 e hp-u), iiij r /2 ^(7)i^or i/2 n < m < i. 

Since X + (7)* = — X_( 7 ), half of (|5.7[) is proved; the other half follows analogously 
from the nonpositivity of Q_. 

For each 7 € D, TV € IN, we define the bounded forms 

00 n in 

q N b)[u,v]~ E ^^( T (l)u,v)\^ , u,veU, 

n=N+l ' ' 

x N (j)[ Ul v] :=q N ("f)[\Ho\ 1/2 u,\Ho\ 1/2 v], u,v e Q. 
Since T( 7 ) is bounded-holomorphic, it clearly holds that 

OO n 

q N ( 1 )[u,v}^(Q N ( 1 )u,v), u,ven, Q N h):= V 1^(0)^0, 

^— ' TV. 

n=N+l 

with convergence in the operator norm topology and uniform in 7 on compact 
subsets of D. Observe that 

f)Lg(7)[^, v] = (P d iag(7K w) - r w (7)[«. «], " e £>(Pdia g ( 7 )), « e Q. 

We now claim that for 7 € (— 1, 1), the operator -ffdiag( 7 ) represents the form 

f)diag(7)M] ~ &(7M7)«,tf(7H «,« e PfoiagW) := U(i)*Q = Q 
and that 2?(7Jdiag(7)) is a core for | _fTo | 1 ^ 2 . This then implies that 
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First, by dill), we have t/( 7 )*Q = Q. Hence, F( 7 ) := \H \ 1/2 U(-/)*\H \~ 1/2 
is onto, one-to-one and bounded and thus an isomorphism by the closed graph 
theorem. It follows that 

|#or /2 Wdia g (7)) = \H \ 1/2 U( 7 )*V(H( 7 )) = n7Mo| 1/2 W7)), 
and the latter is dense in Ji; this proves the second claim. The first claim follows 
from the spectral theorem if we observe that the spectral families {-Ea(7)}a£R and 
{E\(j)}\£R of the self-adjoint operators if (7) and ifdia g (7) are related by 

Exh) = U(^)*E x (j)U(j), A G R. 
In particular, this yields 

(5.9) 2?(|iW 7 )| 1/2 ) = C/( 7 )*P(| J ff(7)| 1/2 ) = £/(7mi#o| 1/2 ) = ^(l#o| 1/2 ), 
whence |if | 1/2 is |fidia g ( 7 )| 1/2 -bounded. Therefore, for u,v e 7i, \\u\\ = \\v\\ = 1, 

|r W (7)[|i?dia g (7)r 1/2 «,| J ffdia g (7)r 1/2 «]| 

=iq^(7)[iHor /2 | J ffdia g ( 7 )r l/2 «j^or /2 ii?dia g (7)r i/2 t']i 

<||Q N (7)|||||i/or /2 |Hdia g (7)|- 1/2 || 2 ^>0 

locally uniformly in 7 e D. Thus, for N sufficiently large, [43J Theorem 2.11] 
yields that there exists a unique self-adjoint operator if <ji ag (7) associated to the 
form f)di ag (7). By construction, T>{H^ ( / y)) is a core for |ffdia g (7)| 1 ^ 2 - Since the 
norms || |ffdia g (7)| 1 ^ 2 • || and || |ifo| 1//2 ■ II are equivalent by (|5.9p . it is also a core 
for I Hoi 172 - 

The norm-resolvent convergence of i?dTa g (7) to ifdia g (7) follows from the final 
argument in the proof of Theorcm l2.6l More precisely, if if C (— 1,1) is a compact 
subset, let 

b N := sup (\\Q N (<y)\\ || |Ho| 1/2 |iW 7 )r 1/2 || 2 ) ■ 

■yeK v ' 

We then have 

IKW7) - ivr 1 (HL e d) - ivr'w < t^-h °- D 

1 - Oat \T}\ 
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